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d _ IRAS 19312+1950 is an SiO maser source that exhibits a prominent bipolar nebulosity. Mapping 

^ '. observations of this object were made in the CO J = 1-0, 13 CO J = 1-0, C 18 J = 1-0, CS J = 2-1, and 

HCN J = 1-0 lines and in the 150 GHz continuum band. Near-infrared imaging observations were also 
made in the J, H, and if-bands. The line profiles of the 12 CO and HCN spectra consist of a weak broad 
component with a line width of about 50 km s _1 and a strong narrow component of the width of about 
3 km s _1 . The profiles of the 13 CO, C 18 0, and CS lines have only the narrow component. Both of the 

■ components have an intensity peak at the IRAS position. The narrow component was clearly resolved 
QO ! with a 15" telescope beam. The spectral energy distribution of this object exhibits a doubly peaked profile 

between 1 and 25 /im. The 150 GHz continuum flux density was found to be 0.07 Jy, which is consistent 
i~ I ' with the flux density predicted by the expanding envelope model with a mass loss rate of ~ 10~ 4 Mq y _1 

at a distance of 2.5 kpc. We argue that the broad component originates from the expanding envelope of 
this object, and that the hot dust cloud, which is the source of the narrow component, is also physically 
associated with this object. Though the present observations do not preclude the possibility of a young 

■ stellar object, we argue that it is less plausible. We conclude that IRAS 19312+1950 is an AGB/post-AGB 
star that is evolved from a massive progenitor. 

q ■ Key words: stars: AGB and post-AGB — stars: circumstcllar matter — stars: individual (IRAS 

19312+1950) 

C/3 '■ 
Ctf ■ 

■ 1. Introduction the central star (~ 10 14 cm), and the OH masers are fur- 
ther away (~ 10 16 cm). Therefore, when the mass loss 

Protoplanctary nebulae (PPNe) arc transient ob- terminates, the SiO masers disappear first, whereas the 
jects between the asymptotic-giant-branch (AGB) and OH masers remain active during the protoplanctary and 
planetary-nebula phases of stellar evolution. They are early planetary- nebula stages (Lewis 1989). Thus, only 
key objects for understanding the formation of planetary a few PPNe with SiO masers have been found in the 
nebulae. Though OH and H2O masers have been ex- past; for example, OH 231.8+1.4 (e.g., Barvainis, Clemens 
tensively searched for in late-AGB stars and PPNe (e.g., 1984; Nyman et al. 1998). Nakashima and Deguchi (2003) 
Likkel 1989; Hu et al. 1994), relatively little work has newly found SiO masers in such cold objects. In particu- 
been made in SiO masers. Observations of OH/IR stars lar, IRAS 19312+1950 has a prominent bipolar nebulosity 
and objects with cold envelopes were made in SiO masers with a large angular size of more than 40" on the 2MASS 
by Jewell et al. (1991), Nyman et al. (1986, 1993), Gomez near-infrared images. 

et al. (1990), and Nakashima and Deguchi (2003). It was The object, IRAS 19312+1950, is unusual in several 

found in these works that the detection rate of SiO and respects. Firstly, it appears to be oxygen-rich (O-rich) 

H2O masers drops with an increase of the IRAS [25-12] according to the detections of SiO, H2O, and OH masers 

color, which is a temperature/evolutionary-stage indica- (Nakashima and Deguchi 2000) and the thermal line of 

tor of the dust envelopes of late- type stars. SO (Deguchi and Nakashima 2003). But, in addition, 

Though H2O and OH masers are often found in PPNe, the thermal line of H 13 CN was also detected toward this 

SiO masers are rarely detected in these objects. The SiO object (Nakashima and Deguchi 2000). This molecule is 

masers are emitted in circumstcllar shells very close to normally considered to be an indicator of carbon-rich (C- 

Z Z 7 j 7^ I r « , TT . .. rT11 . . rich) circumstcllar chemistry under the equilibrium con- 
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61801-3080, U. S. A. non-equilibrium (Deguchi, Goldsmith 1985; Deguchi et 
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al. 1986). Secondly, a doubly peaked profile was observed 
in the spectrum of the SiO maser in the J =1-0, v ~2 
transition, and SiO maser emission exhibits a strong time 
variation (Deguchi and Nakashima 2003). The spectral 
features of OH masers are also quite unusual, exhibiting 
multiple spiky peaks around the systemic velocity (B. M. 
Lewis 2000, private communication). 

This object has a cool dusty envelope; the IRAS color, 
log(F 2 5/Fi 2 ) = 0.47, indicates the presence of cool dust 
with a temperature of about 200 K around the central 
star. This IRAS color is at the boundary between AGB 
and post- AGB phases in the two-color diagram (Nyman ct 
al. 1998). Though we postulated that this is a protoplan- 
etary nebula, it is reasonable to think if this is a young 
stellar object with SiO masers in a molecular cloud, like 
IRc 2 in Orion. In fact, Nakashima and Deguchi (2000) 
discussed a double-peak SiO maser profile of this object, 
which is relatively rare in AGB stars, but is well known 
in Ori IRc 2. However, based on the various observational 
facts that have been collected since the discovery, we be- 
lieve that this is an evolved object. A discussion on this 
point will be given in section 3 of this paper and in a 
forthcoming paper. 

Because IRAS 19312+1950 has both prominent and un- 
paralleled characteristics as a protoplanetary nebula, we 
decided to observe this object more carefully. This object 
may provides an important clue to solve a puzzle of mass 
ejection mechanism in a class of cool objects. In this pa- 
per, we report on the result of mapping observations in 
molecular lines of 12 CO, 13 CO, C 18 0, HCN, and CS, and 
in the 150 GHz continuum band using the Nobeyama 45- 
m telescope. The spatial and kinematic structures of the 
dust /molecular envelope around the central star have been 
revealed. In addition, we made photometric and imaging 
observations in the near-infrared J, H and K bands. We 
also made molecular-line observations by the other lines, 
which clarify the special circumstellar chemistry of this 
object. A part of the circumstellar chemistry of this ob- 
ject will, however, be presented in another paper. 

2. Observations and Results 

2.1. Observations 

Mapping observations in the CO J = 1-0, 13 CO J = 1 
0, C 18 J = 1-0, CS J = 2-1, and HCN J = 1-0 lines 
were made with the 45-m telescope at Nobeyama during 
2001 February 4-9. A cooled SIS-mixer receiver with a 
bandwidth of about 0.5 GHz was used, and the system 
temperature (including atmospheric noise) was 300 - 500 
K (SSB), depending on the frequency and weather. The 
aperture efficiencies at 86, 100, and 110 GHz were 0.41, 
0.36, and 0.31, respectively. The main beam efficiencies 
at 86, 100, and 110 GHz were 0.49, 0.45, and 0.42, respec- 
tively. The half-power full widths of the main beam were 
19". 8, 17". 1, and 15".6 at 86, 100, and 110 GHz, respec- 
tively. The antenna temperature (Ta) given in the present 
paper was corrected for the atmospheric and telescope 
ohmic loss, but not for the beam or aperture efficiency. 
Acousto-optical spectrometers with wide frequency cover- 



age (AOS-W) were used. The AOS-W spectrometers have 
250 MHz bandwidth and 2048 effective channels, the ve- 
locity coverage of about 800 km s _1 and a spectral resolu- 
tion of about 0.85 km s _1 at 86 GHz. The rest frequencies 
of observing molecular lines were taken from Lovas (1992). 
For an intensity calibration and a check of tuning frequen- 
cies, we observed molecular lines toward W51 (e.g., Bell et 
al. 1993) prior to the observations of IRAS 19312+1950, 
and confirmed that the strong lines were surely detected 
toward this cloud. The telescope pointing was checked 
using a strong SiO maser source, Villi Oph, and it was 
confirmed to be better than 3 under windless conditions. 
Mapping patterns for each line are summarized in table 
1. The bipolar axis of this nebula seen on the 2MASS 
images is oriented approximately by 45° from the north 
(to the east). Therefore, cross and rectangular mappings 
were made along the perpendicular and parallel directions 
to the bipolar axis. The IRAS position of this object (19 
31 m 12 s . 8, 19° 50' 22", B1950) was set as the map centers. 
Since contamination by background molecular emission 
was expected, we carefully chose the offset position for po- 
sition switching; it was 20' away from the mapping center 
(in the galactic latitude). The raw data were processed, 
flagging out bad scans, integrating with r.m.s. -weighted 
means, and removing the slopes in baselines. 

Table 1. Mapping patterns 
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The continuum observations at 150 GHz were made on 
2001 January 11 and 12, with the Nobeyama Bolometer 
Array (NOBA) mounted on the Nobeyama 45-m tele- 
scope. The NOBA is a seven-element bolometer array 
using a differential readout method (Kuno et al. 1993). 
The band pass of the filter was centered at 147 GHz with a 
bandwidth of 30 GHz. The half-power full width (HPFW) 
of the beam was 12". Maps were made with raster scans 
along the right ascension and declination. The scanning 
speed was 15" s . These maps were combined using the 
basket-weaving method to correct scanning effects. The 
angular size of the final map is 1'.9 x 1'.9 with a grid spac- 
ing of 5". 3. The final integration time per sampling point 
was about 40 s. The r.m.s. noise level of the final map is 
~ 10 mJy beam -1 . Further details of the NOBA are given 
in Kuno et al. (1993). Pointing and flux calibrations were 
performed every hour by observing 3C 345 and NGC 7027. 
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Table 2. Observed line intensities for IRAS 19312+1950 
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T & y. peak intensity of original spectrum. 

T & 2'- peak intensity of background-subtracted spectrum. 

J.J.: integrated intensity of background-subtracted spectrum. 



We assumed that the spectrum of 3C 345 is flat in the ob- 
serving bandwidth. The effect of atmospheric extinction 
was corrected using the opacity of the atmosphere mea- 
sured between observations. According to the uncertainty 
of these calibrations and the variation in the antenna ef- 
ficiency due to wind, the uncertainty of the absolute flux 
density was estimated to be ±30%. 

Near-infrared photometric observations were made with 
the 2.3-m telescope of Australian National University at 
the Siding Spring Observatory, Australia, on 2000 July 7, 
using the infrared array camera CASPIR. The detector of 
the CASPIR was a 256x256 InSb array with the angular 
size of 0".5 per pixel at the sky, and it covers a field of view 
of 128" x 128". The seeing size was l"-2". In order to pre- 
vent the saturation of the array, a short exposure time was 
used (0.3 s x 100 cycles for the H and K bands, 1 s x 100 
cycles for the J band, and 12 s x 10 cycles for the B 7 and 
H2 narrow bands). To determine the photometric magni- 
tudes in the standard system, we took images of standard 
stars near IRAS 19312+1950. The data were reduced us- 
ing the IRAF data-reduction package by a standard way: 
dark frames were subtracted from raw data; flat fielding 
was done with dome-flat frames and the magnitudes were 
measured by an aperture photometry technique using the 
DIGIPHOTO package in the IRAF. The raw images us- 
ing the H2 narrow-band filter (2.12 fim) exhibited a faint 
patchy structure a few second north of the central star. 
However, in the later analyses, this was turned out to be 
a ghost of the bright central star produced by reflection 
in the narrow-band filter. We detected no H2 emission in 
this nebula. 

2.2. Results of Molecular Line Mappings 

Figure 1 shows the spectra of the CO J = 1-0, 13 CO 
J = 1-0, C 18 J = 1-0, HCN J = 1-0, and CS J = 2-1 
transitions toward IRAS 19312+1950. The line parame- 
ters are given in table 2. We detected strong CO emis- 
sions toward IRAS 19312+1950. However, the 60"-grid 
9-point cross mapping revealed that the emissions peaked 
at V| sr =27-33 and 43 km s _1 are spatially extended more 
than a few arc minutes, indicating that these are from 
back/foreground molecular clouds. With the 15" grid 
mapping, we found that the component at Vj sr = 38 km 
s _1 is sharply peaked in intensity at the IRAS position. 
The raw spectra indicate that the spectral features of 



back/foreground clouds are relatively homogeneous over 
the mapping regions. Therefore, in order to obtain the 
spectra of the object, we averaged the spectra at the 22" 
to 60" offset positions surrounding the object, and made 
an average background spectrum. We then subtracted 
it from the original on-source spectra and obtained the 
(supposedly) background-free spectra. The same proce- 
dures were applied to all of the lines, though contami- 
nation by background emission was not so severe except 
in the cases of the CO J = 1-0 and 13 CO J = 1-0. The 
fore/background components of the CO J = 1-0 line were 
considerably flat and the intensity varies by only less than 
10% in the 1' away spectra. In the case of 13 CO J = 1-0, 
the background emission was not perfectly flat, so that the 
subtracted spectrum still exhibits emissions at V[ gr =28 
and 45 km s _1 . Figure 1 shows the raw spectra at the 
mapping center (thin line), the average background spec- 
tra (broken line), and the background-subtracted spectra 
at the mapping center (bold line). In further analyses, 
we used the background-subtracted spectra to avoid the 
background contamination. 

All of the spectra exhibit a clear narrow feature with 
a width of about 2-3 km s~ x (hereafter narrow compo- 
nent) around V[ sr = 36 km s _1 . The widths of the narrow 
components seem to vary among the molecular lines. For 
instance, the lines of 12 CO, 13 CO, and C ls O have full 
widths at half maximum of about 3.0, 2.3, and 2.0 km 
s" 1 , respectively. 

On the other hand, the background-subtracted spec- 
trum of 12 CO (solid line in top left panel of figure 1) ex- 
hibits a weak broad feature with a line width of about 
40 km s" 1 . This was quite weak and difficult to distin- 
guish from noise or background emission. However, in the 
top- left and right panels of figure 1, we can clearly see 
emission at Vj sr = 20-25 and 48-55 km s _1 , which are sig- 
nificantly above the background emission (broken curve) 1 . 
The presence of this component can be seen more clearly 
in the HCO + spectrum, which is free from background 
contamination (Deguchi and Nakashima 2003). The broad 
component was also seen in the spectra of thermal lines 
of SO and SiO (Deguchi and Nakashima 2003). In the 

1 Though the presence of the broad component seems to remain 
somewhat debatable in this single-dish observation, we have con- 
firmed it with mm-wave interferometric observations. The re- 
sults described in this paper are completely consistent with the 
interferometric observations. 
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case of HCN, we can also see a possible broad component, 
though it is not so clear because of the hyperfinc splitting 
of the HCN J = 1-0 transition: F = 0-1, 2-1, and 1-1. 
The hyper fine lines should appear at AV = —7.44 and 
+4.68 km s" 1 relative to the main F = 2-1 line (89.08792 
GHz) with an LTE (optically thin) line intensity ratios 
of 1:5:3 (e.g., Cao et al. 1993). In the bottom left panel 
of figure 1, the origin of the radial velocity of the HCN 
spectrum is based on the rest frequency of the strongest 
transition F = 2-1. 

Figure 2 shows the integratcd-intensity maps. In this 
figure, the F-axis is taken along the bipolar axis of the 
nebula; it is inclined by 45° from the north. The back- 
ground emission was subtracted in these maps. For mak- 
ing these maps, the data taken on different days (with 
different mapping patterns) were combined, because the 
wind often prevented us to continue mappings (of course, 
the data taken under the wind speed more than 10 m s" 1 
were discarded). In the cases of the 12 CO and HCN lines, 
we made two separate maps: one for the broad and an- 
other for the narrow components. We took the emission in 
the radial velocity ranges from to 30 km s _1 and from 
50 to 60 km s _1 for the broad component of the 12 CO 
spectrum, and from to 25 km s , and from 45 to 60 
km s _1 for the broad component of the HCN. The veloc- 
ity ranges of the broad component in the HCN spectrum 
were chosen by eye inspection. The other lines (except for 
12 CO and HCN) do not exhibit the broad feature; maps 
were made only for the narrow components. 

A remarkable point in figure 2 is that all of the intensity 
peaks of the narrow and broad components of molecules 
coincide with the IRAS position of IRAS 19312+1950. We 
do not find any differences in the peak position between 
the narrow and broad components, or between different 
molecular lines. The object seems to be slightly resolved 
by the 15" beam. The directions of the elongation seem to 
be different between molecular lines; the narrow compo- 
nents of the 12 CO and HCN spectra are somewhat elon- 
gated along the X-direction; the 13 CO and C 18 spectra 
are elongate along the ^-direction; the broad component 
of the 12 CO spectrum. 

In order to evaluate the size of the source more carefully, 
we computed the model intensity variations and compared 
these with the observations. Figures 3a and 3b show po- 
sition vs. integrated-intensity plots along the X and Y 
directions for each line. The intensities expected from 
simple homogeneous circular-disk models with radii of 0" 
(point source), 10", 20", and 30" are drawn in solid and 
broken curves alternatively The horizontal and vertical 
tick marks indicate the uncertainty in telescope pointing 
and the rms of integrated intensities, respectively. In fig- 
ure 3a, we can see that the extension of the narrow com- 
ponent of 12 CO is clearly larger than that of the broad 
component. However, this phenomenon is not so clear in 
the case of HCN. It seems that the narrow components of 
all of the lines, except for that of HCN, extend to more 
than 15" according to the simple homogeneous disk model. 

Figures 4a and 4b show position-velocity (PV) dia- 
grams for each molecular line (background subtracted). 



No particular spatial velocity gradient was seen in the 
12 CO and CS. The narrow components (~ 35-36 km s _1 ) 
of 13 CO, C 18 0, and HCN exhibit a small velocity gradi- 
ent along the X-axis (3.5 x 10 -2 km s _1 arcsec -1 ) that 
is perpendicular to the bipolar axis. Even though the 
background emission was not completely homogeneous, 
the intensities of the background clouds at around 30 and 
45 km s _1 remained to be small compared with the nar- 
row components, especially in the 12 CO, 13 CO, and C 18 
lines. 

2.3. Results of 150 GHz Continuum and Near- Infrared 
Observations 

Figure 5 shows the 150 GHz continuum map toward 
IRAS 19312+1950. The map was overlaid on the JHK 
false-color composite image. The image size is 1'.75 x 
1'.75. The 150 GHz continuum emission seems to extend 
to the north-east along the bipolar axis. The peak flux 
density of the 150 GHz continuum radiation was 7.7 x 10 -3 
Jy per beam. The total integrated flux density of the 
source integrated over the inner 25" was 0.07 Jy (±0.01 
Jy) . The near- infrared image shows that the number den- 
sity of bright stars is less at the northeast part of the 
object than at the southwest part, suggesting that the 
dust hides background stars more in the northeast part of 
the object. The magnitudes of the central star (excluding 
the ring-like elongated component) are 10.7, 7.8 and 6.3 
mag (uncertainty: ±0.1 mag) in the J, H, and iC-bands, 
respectively. Figure 6 shows the spectral energy distri- 
bution (hereafter SED) of IRAS 19312+1950. Black body 
curves are also drawn. The flux densities of 12, 25 60, and 
100 /im were taken from the IRAS point source catalog. 
This SED shows a flat (or moderate double peak) shape. 
An excess can be seen between 60 and 100 /im. These 
two data points were fit by the blackbody curve of the 
temperature of 70 K. In the galactic plane, interstellar 
dust in background molecular clouds occasionally emits 
far-infrared radiation. Therefore, this excess may possi- 
bly be caused by interstellar contamination. However, the 
IRAS mid- and far-infrared maps with about 1 resolution 
show a clear intensity peak at the IRAS position of IRAS 
19312+1950. This fact suggests that most of the energy 
of radiation in the 60 and 100 /im bands originates from 
the central source of IRAS 19312+1950. 

3. Discussion 

3.1. Evidence of an Evolved Object 

Because molecular species observed toward this source 
are somewhat similar to the molecules detected in dark 
clouds (for example, in L 134N; Swade and Schloerb 
1992), it is reasonable to doubt if this star is a young stel- 
lar object with SiO masers as IRc 2 in Orion. However, we 
think that this is not the case for IRAS 19312+1950. Note 
that many molecular lines have also been detected in the 
protoplanetary nebula, OH 231.8+4.2 (Morris et al. 1987). 
Furthermore, AGB maser stars have occasionally been 
found in/toward molecular clouds (Imai et al. 2002) 

Firstly, this object exhibits deep 2.3 /im CO first- 
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overtone absorption bands corresponding to late M III 
(P. R. Wood 2002, private communication). Most of low- 
mass young stellar objects exhibit emission or weak ab- 
sorption of the 2.3 fim CO overtone bands; this fact is 
very often used to separate YSOs and AGB stars (Casali 
and Matthews 1992; Chandler et al. 1993). Exceptional 
cases are IRc 2 in Orion and IRS 5 in L 1551, which ex- 
hibit the 2.3 /jm CO absorption spectra somewhat similar 
to those of K I and K III, respectively (Morino et al. 
1998; Carr et al. 1987). However, even in these excep- 
tional YSOs, the effective temperatures are much higher 
than those of M III stars. 

Secondly, IRAS 19312+1950 is an isolated object. We 
find no bright object within a few arcminutes from this 
star in the optical, near-infrared, and middle-infrared 
(MSX) images. Moreover, we cannot find any star form- 
ing activities on the infrared images down to K ~ 16 mag 
(note that interstellar extinction at K band in this region 
is about 2 mag at most, which is estimated in subsec- 
tion 3.2). The closest IRAS source (IRAS 19309+1947) is 
located 5' south-west. We find several very faint MSX ob- 
jects at about 4' - 10' away and radio continuum sources 
(probable H II region) about 10' north-west (with the 
SIMBAD database). These objects are located too far (> 
3 - 7 pc at the distance of 2.5 kpc) to be considered as 
a cluster of star-forming activities associated with IRAS 
19312+1950. They are probably forc/background objects. 

Because both of above two reasonings are still not com- 
pletely crucial as evidence of an evolved star, we do not 
preclude the possibility of a young stellar object here. 
However, because of the absorption feature of CO over- 
tone bands, the central star of IRAS 19312+1950 cannot 
be a low-luminosity young stellar object. If we consider a 
high-luminosity 10 4 L Q ) young stellar object, the dis- 
tance and physical conditions must be somewhat similar 
to those of AGB/post-AGB stars. Therefore, we postu- 
late that it is an AGB/post-AGB object in this paper, 
and estimate various quantities. Further discussion on 
the possibility of a young stellar object will be given in 
a future paper, which is based on the molecular species 
found in this object. 

3.2. Distance to the Object 

We can estimate the distance to the source from the 
SED of this object, applying an interstellar extinction 
law at each wavelength (for example, see Mathis 1990). 
However, the distance depends strongly on the amount 
of total extinction toward this particular object. The ex- 
tinction of the CO fore/background cloud at V| sr = 27 - 
33 km s _1 is estimated to be Ay = 16.3 from the CO 
integrated intensity (24.7 km s ; for conversion factor, 
see Frcrking et al. 1982). If this value is used for the 
extinction assuming that this is a foreground cloud, we 
obtain the distance of 2.5 kpc from the absolute luminos- 
ity of the object, 8000 L@, which is typical for AGB stars. 
The bolometric correction to the IRAS 12 /im flux density 
(Deguchi et al. 1998) was computed to be (BC)i2 = 3.8, 
where halves of the IRAS flux densities at 60 and 100 fim 
are assumed to come from this object. 



If we consider that the V[ sr = 27-33 km s" 1 component 
cloud is behind the object (though this is unlikely) and 
apply the standard interstellar extinction (for example, 
see Deguchi et al. 1998), we obtain the distance of 3.9 
kpc. The distance varies from 5.1 to 3.3 kpc, when all of 
the 60 and 100 /xm flux densities come from the source, or 
nothing comes, respectively. 

The luminosity of the central star of PPNe is known 
to be larger than the assumed value, 8000 L Q . In such a 
case, the distance varies in proportion to the square root 
of the luminosity. On the other hand, the distance would 
be slightly smaller than the above value if we take into 
account the self-absorption due to dust in the narrow- 
component cloud. Therefore, we estimated that the dis- 
tance ranges from 2.0 kpc to 5.1 kpc for this object. If 
we compare this with OH 231.8+4.2 (at the distance of 
1.3 kpc; Bowers, Morris 1984), IRAS 19312+1950 is still 
a bit a distant object, though the IRAS flux densities at 
middle infrared wavelengths arc similar for both objects. 

For simplicity, we assume in this paper that the distance 
to the object is 2.5 kpc, and estimate various quantities. 

3.3. Broad Component in the 12 CO Spectrum 

In the background-subtracted spectrum of 12 CO (top 
two panels of figure 1), a broad component is seen in the 
velocity range V\ sr ~ 10 - 60 km s _1 . This component 
is considered to be formed in the expanding envelope of 
the central mass-losing star. Because the broad compo- 
nent is weak and overlapped with the strong background 
emission, it is difficult to measure the accurate full line 
width of this weak component. The full width at zero 
maximum is estimated to be about 50 km s , suggest- 
ing a typical expanding velocity of the envelope of 25 km 
s . In addition, the velocity range of SiO maser emis- 
sion, which is found in Vi sr =20-55 km s _1 (Nakashima 
and Deguchi 2000), is overlapped with the radial veloc- 
ity of the CO broad component. Moreover, the intensity 
map of the broad component (see top left panel of figure 
2) shows a clear intensity peak at the center position of 
this object. These facts strongly suggest that the broad 
component is physically associated with the central star. 
From the CO line intensity (and the intensities of other 
molecular lines), the mass loss rate can be calculated with 
the LVG model (for example, see Morris and Alcock 1977). 
The best-fit value of the mass loss rate was calculated to 
be 1.04 x 1O~ 4 M yr -1 . The details of the LVG model 
calculation are to be given in a future paper. 

3.4- Narrow Components 

We should be careful in interpreting the narrow com- 
ponent. The width of the narrow component is about 
3 km s _1 ; this is nearly the same for all of the lines 
of 12 CO, 13 CO, HCN, and CS. On the other hand the 
widths of the narrow components of less abundant species 
(such as HCO + , etc.) are dependent on molecular species 
(Deguchi and Nakashima 2003). This fact suggests that 
these molecular lines arc formed in a turbulent molecu- 
lar cloud (e.g., Lada 1976); the line width depends on 
the optical depth of the line in such a cloud. The nar- 
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row component emission may come from some cool com- 
ponent toward this object, for example, a dark cloud at 
fore/background of this source. 

We can estimate a column density of H2 from the CO 
line intensity when we assume that the narrow component 
of CO comes from a dark cloud. Adopting the CO-H 2 
conversion ratio of 2.6 x 10 20 /(K km s -1 ) (Magnani et 
al. 1998), we obtained a H2 column density of 6.1 x 10 21 
cm~ 2 from the present observation [Tp ea j c ( 12 CO)=3.9 K, 

HPFW=2.9 km s _1 , and the telescope beam efficiency 
of 0.48]. The total mass of H2 within the 16" telescope 
beam (assuming the distance of 2.5 kpc) is calculated to be 

2.8 Mq. Alternatively, from the column-density estimate 
using optically thin C 18 line (for example, see Frcrking 
et al. 1982), we obtain the mass of about 9.7 Mq for the 
cloud within the 16" telescope beam. The mass estimates 
using CO lines arc considered to be uncertain by a factor 
of few, depending on the isotopic abundance ratio. 

In the P-V diagram (figures 4a and 4b), the narrow 
components of 13 CO and C 18 exhibit a small velocity 
gradient along the direction perpendicular to the bipo- 
lar axis. We estimated the mass of the source based on 
this velocity gradient, assuming the Keplerian rotation of 
the narrow component, 0.5 km s _1 at an inner radius of 

1.9 x 10 17 cm (5"). Then, we obtained about 3.5 M as 
a mass of the object inside the radius of 7". In fact, an 
inclination of the ring structure [which is seen with the 
SUBARU and UH88 telescopes (M. Tamura 2001, private 
communication)] is estimated to be ~ 67°. In this case 
the mass would increase by a factor of 1.2. 

3.5. Dust Emission 

We detected continuum emission at a level of 0.07 Jy 
at 150 GHz. This emission is considered to be radiation 
emitted by dust grains toward this object. 

We first estimate the mass of the cloud from a simple 
homogeneous-cloud model (e.g., Hogerheijdc et al. 1997). 
When we assume = 10 K and a distance of 2.5 kpc, 

we obtain the optical depth of the dust to be 1.4 x 10 -3 
at 150 GHz within the telescope beam. Adopting a mass- 
absorption coefficient (k„) of 3 x 10 -3 cm 2 g _1 (including 
both gas and dust particles), we obtain 140 Mq for the 
total mass of the dust cloud within the telescope beam. 
This value is inconsistent with the mass estimated from 
the CO narrow component, suggesting that the 150 GHz 
continuum emission docs not come from the cold dark 
cloud. 

Because the dust cmissivity increases with the temper- 
ature of the grain, it is possible that the dust temperature 
is much higher than 10 K. In fact, the far infrared spec- 
trum in figure 6 can be fit by a blackbody with Td us t=70 
K or higher. The dust temperature of an expanding enve- 
lope of a mass-losing star is higher near to the central star 
and decreases with the radius. Knapp et al. (1993) cal- 
culated the continuum flux densities of stars with various 
mass-loss rates, and found that the model flux densities 
fit well with the observed spectra. Their model for OH 
231.8+1.4 with a mass-loss rate 2 x 10~ 4 Mq yr _1 gave 



^150 GHz = 0-39 Jy at a distance of 1.3 kpc. If the same 
model is applied to IRAS 19312+1950, the flux density 
at 150 GHz is to be 0.1 Jy at the distance of 2.5 kpc. 
This value is roughly consistent with the continuum flux 
density that was obtained in the present paper. In this 
model, the total mass of the gas in the envelope within 
the radius of 9.3 x 10 17 cm (25") is computed to be about 
0.1 Mq (corresponding to the broad component). 

The dust is also responsible for the nebulosity seen on 
the near-infrared images, which is made by scattering the 
light from the central star. We can estimate the amount 
of dust grains by assuming that the dust optical depth at 
the K band is about unity in the envelope of this star. 
In this case, a simple estimate of the mass of the the ring 
structure, which is present at the radius of 5", gives 0.9 
Mq. Here, we assumed a l"-thickness spherical dust shell 
and the standard gas-dust ratio. 

From these estimates of the masses in this and previ- 
ous sections, and from the fact that the observed velocity 
of the narrow component (~35-37 km s" 1 ) coincides well 
with the center velocity of the broad components, we con- 
clude that the cloud responsible for the narrow compo- 
nent is not physically independent from the central star. 
Rather, we are compelled to consider that the narrow com- 
ponent cloud of a few solar mass is a part of the object. 

3.6. Is IRAS 19312+1950 a Protoplanetary Nebula ? 

In general, PPNe are expected to have the following 
observational characteristics (e.g., Kwok 2000); (1) strong 
infrared excess that is due to circumstellar matter; (2) 
the effective temperature of the central star between F 
and G types, (3) the presence of a detached dust shell 
due to the termination of an AGB-phase mass loss, which 
gives a doubly peaked spectral energy distribution; (4) the 
nebulosity seen in scattered light and no line emission of 
highly ionized atoms. 

IRAS 19312+1950 shares the first and the forth proper- 
ties. It exhibits a red IRAS color [log(F 25 /Fi 2 ) = 0.47] and 
shows a clear intensity peak at the IRAS position in the 
12 CO intensity map. The B 7 narrow-band imaging with 
the ANU 2.3-m telescope resulted no evidence of ionized 
gas. At mid-infrared wavelengths, no line emission fea- 
ture was recognized in the IRAS low-resolution spectrum 
(except a probable noise at 17.8 /im; Volk et al. 1991). 
On the other hand, we should be careful in characteris- 
tics (2) and (3). Effective temperatures of stars between 
AGB and PN phases are 5000 - 7000 K, which correspond 
with the surface temperatures of F- or G-type stars. In 
many cases, the central stars can be observable at optical 
wavelengths. A quite faint counterpart is recognizable at 
the position of the central star in the phase II DSS image 
in R-band, but the star image is too faint to judge the 
spectral type. One explanation of this discrepancy in the 
characteristics is that IRAS 19312+1950 suffers a large 
interstellar extinction of Ay ^16 as estimated in subsec- 
tion 3.2. The spectral energy distribution does not show 
clear single nor double peaks. However, if the interstellar 
extinction (Ak = 1.83) is corrected, the fluxes at 0.9 - 
2.2 /j,m are increased by a factor of 30 - 7, respectively. 
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Therefore, the spectrum seems to be doubly peaked. The 
J-, H-, and K- band intensities are well fit by a black- 
body with T = 1770 K when the extinction of A K = 1.83 
is corrected. The flux densities at 12 and 25 /zm are well 
fit by a blackbody with T = 190 K, which is a temperature 
reasonable for the dust shell at medium radii in the enve- 
lope (Hashimoto 1994). These properties also suggest that 
IRAS 19312+1950 is a very young protoplanetary nebula. 

Progenitors of O-rich protoplanetary nebulae are in- 
ferred to be relatively massive stars (> 4 Mq). According 
to the evolutionary scenario at the AGB phase, the stars 
with masses below ~ 4 Mq ends up as a C-rich star 
by the He-burning at the core and succeeding dredge-up. 
However, in stars above that, the created 12 C atoms are 
further fused by the Hot Bottom Burning (HBB: Iben, 
Renzini 1983; Wood et al. 1983) at the late-AGB phase 
and the massive stars remain to be O-rich throughout the 
AGB phase with experiencing no (or a very short) C-rich 
era (Marigo 2001). Therefore, the mass of the progenitor 
of IRAS 19312+1950 (O-rich) is estimated to be relatively 
large. Recently, Ventura et al. (1999) suggested that stars 
with a relatively large initial masses (~ 3.8 Mq) first go 
through a short (At ~ 30000 yr) J-type ( 13 C-enhanced) 
carbon star era before the hot bottom burning converts 
them back into oxygen-rich stars, and ends back again 
as C-rich stars. The estimated mass of the central star 
of IRAS 19312+1950 is somewhat too large to experience 
this evolutionary scenario. 

The coincidence of the masses, about a few Mq, which 
were calculated for the dust shell and the CO narrow com- 
ponent in the previous sections, even suggests us another 
idea that IRAS 19312+1950 is the AGB object embedded 
in a small globule-like cloud. The near-infrared radia- 
tion emitted by the central star can illuminate the dusty 
cloud having a mass of about a few Mq. In this case, 
the cloud cannot be a cradle of the AGB progenitor be- 
cause the mass of the progenitor is comparable with this 
cloud. Only one possible case would be the AGB star col- 
liding with the cloud (for example, a cometary cloud; see 
Kulkarni et al. 1992]. Even in such situation, the relative 
velocity between the cloud (the narrow component) and 
the AGB star (the broad component) is so small (< 1 km 
s _1 ) that it is unlikely for such event to occur often. 

The line profiles of CO with broad and narrow fea- 
tures have been found in many stars (Knapp et al. 1998; 
Kerschbaum and Olofsson 1999). These features are in- 
terpretated by multiple winds or Keplerian disk (Bergman 
et al. 2000). Therefore, it is possible to consider that the 
narrow-component cloud in IRAS 19312+1950 is mate- 
rial trapped in an orbiting molecular reservoir (Jura and 
Kahanc 1999). It has been suggested that the origin of 
bipolar nebulae is a binary system [Morris 1981; Cohen 
1981; Sanchez Contreras et al. 2000]. However, in this 
case, it also seems difficult to keep the ejected material of 
a few solar masses moving around the central star; precise 
modeling of this object is necessary. Because the trapped 
material (supposedly the ring structure) is located consid- 
erably far (5") from the central star, the secondary star 
must be located also at the similar distance. However, 



we find no evidence for the other late-type star at such 
distances. Rather, it is possible that the counter part 
is evolved quickly to a white dwarf, and remnant of the 
ejected material still circles around the binary system, as 
suggested for the silicate carbon stars, etc. (Lloyd Evans 
1990; Jura et al. 1995; Kahane et al. 1998). The exis- 
tence of a hot secondary companion, however, enables us 
to make models explaining the complex chemistry of IRAS 
19312+1950. 

With these arguments, we think at present that IRAS 
19312+1950 is likely to be at an early stage of protoplan- 
etary nebula. However, the origin of the narrow com- 
ponent is still a puzzle. The chemistry of this compo- 
nent is quite similar to that of dark clouds (Deguchi and 
Nakashima 2003) , rather than that of protoplanetary neb- 
ula (Morris et al. 1987), though the dust temperature 
seems too high (^70 K). They may be a remnant of ma- 
terial, that is ejected in the AGB phase from the cen- 
tral star and is trapped by a distant binary counterpart. 
Alternatively, a more drastic hypothesis may be a rem- 
nant of star forming material left after an encounter with 
a molecular cloud. More discussions on these points will 
be given in the future. 

4. Conclusion 

We made the mapping observations of IRAS 
19312+1950 in molecular lines of CO J = 1-0, 13 CO 
J = 1-0, C 18 J = 1-0, CS J = 2-1, and HCN J = 1-0 
lines and in the 150 GHz continuum. We also made 
photometric observations in the J, H, and K bands 
to obtain the spectral energy distribution of IRAS 
19312+1950. The main results are as follows: 

• The line profile of 12 CO (and possibly HCN) exhibits 
weak broad (~ 40 km s _1 ) and strong narrow (~ 
3 km s _1 ) components, whereas the profiles of the 
other lines ( 13 CO J = 1-0, C ls O J = 1-0, CS J = 2 
1) exhibit only narrow components. Both of the 
components show a clear intensity peak at the IRAS 
position of IRAS 19312+1950. 

• The narrow component spatially extends at least 
more than 15 . 

• The spectral energy distribution of IRAS 
19312+1950 exhibits a moderate double-peak 
profile (if extinction is corrected) at the near- and 
mid-infrared wavelengths with excess at 60 and 100 
(im. 

Based on the present observational results, we suggested 
that IRAS 19312+1950 is a transient object between the 
AGB and planetary nebula phases, possibly a protoplane- 
tary nebula. Though the present observations do not com- 
pletely exclude the possibility that this object is a young 
stellar object, evidence of an evolved object is accumu- 
lating furthermore with mm-wave interferometric obser- 
vations. 

The authors thank staff of Nobcyama radio observatory 
for the help of observations. They also acknowledge Peter 
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= 1-0 (top left), enlargement of 12 CO J = 1-0 (top right), 13 CO J = 1-0 (middle left), C ls O J = 1-0 
—0 (bottom left) and CS J = 2-1 (bottom right). The bold line indicates a background-subtracted 
indicates raw spectrum at the source center. In the panels of CO and 13 CO, the broken line indicates 
the spectra at (X,Y) = (0", ±60"), (±30", ±30"), (0", ±30"), (±60", 0") aid (±30", 0"). The spectra 
±60"), (±60", 0"), (±30", 0") and (0", ±30") for ls CO, at (0", ±60"), (±1 5", ±15") and (±60", 0") 



Fig. 1. Spectra of i2 CO J 
(middle right), HCN J 
spectrum, and the thin lin< 
the spectrum averaged over 
were averaged over at (0" 
for HCN, and at (0", ±60"), (±30", ±30"), (0", ±30"), 
indicates the best fit parabola to the broad component. 



(±60", 0") and (±30", 0") for CS. The broken line in the top right panel 
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Fig. 2. Integrated intensity maps of the 12 CO, 13 CO, C ls O, HCN and CS molecular lines. The broad (left) and narrow (right) 
components of 12 CO and HCN are indicated by (B) and (N). Observed positions are indicated by the small dots. The Y-axis is 
taken along the bipolar axis of the nebula, i.e., in the direction 45° from the north, and the X-axis that is perpendicular to it. The 
contour interval and the peak intensity are 1.0 and 14.3 K km s — 1 , 1.0 and 18.3 K km s — x , 2.0 and 29.2 K km s , 0.5 and 5.7 K 
km s -1 , 0.5 and 5.9 K km s -1 , 0.1 and 3.2 K km s" 1 and 0.5 and 8.6 K km s" 1 for CO (B), CO (N), 13 CO, C ls O, HCN (B), HCN 
(N) and CS, respectively. The beam sizes of the 45m telescope at 89 and 115 GHz are given in the bottom right panel. 
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Fig. 3a. Spatial intensity variation of the broad and narrow components in the 12 CO (left) and HCN (right) spectra. The filled 
and open circles indicate the data at offset positions along the X and Y directions, respectively. The solid and broken curves show 
the intensity variations that are expected from the simple circular disk model with radii of 0" (point source), 10", 20" and 30", 
alternatively. In the model, we assumed that the beam pattern of the telescope is Gaussian. The peak values of the model curves 
are adjusted to the observed values. 
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Fig. 3b. the same as figure 3a but the narrow components of 13 CO (top), C ls O (center) and CS (bottom). 
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Fig. 4a. Position-velocity diagrams for the 12 CO (upper two panels), 13 CO (Center two panels) and 18 CO J = 1—0 (lower two 
panels) spectra along the X (left) and Y (right) axes. The contour intervals and the peak intensities are 0.4 K and 3.4 K for 12 CO, 
0.5 K, and 3.4 K for 13 CO and 0.1 K and 0.9 K for C ls O. 
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Fig. 4b. Same diagrams with figure 4a but for HCN (upper two panels) and CS (lower two panels). The contour intervals and the 
peak intensities are 0.05 K and 0.5 K for HCN and 0.1 K and 0.9 K for CS. 
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Fig. 5. 150 GHz continuum map (contour) toward IRAS 19312+1950 overlaid with the composite color image of the J, H and K 
bands (epoch: B1950). 
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